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We study spectral statistics of a Gaussian unitary critical ensemble of almost di- 
agonal Hermitian random matrices with off-diagonal entries ~ — 
small compared to diagonal ones ~ 1. Using the recently suggested method 
of virial expansion in the number of interacting energy levels (J.Phys.A 36,8265 
(2003)), we calculate a coefficient oc 6^ ^ 1 in the level compressibility x{b)- We 
demonstrate that only the leading terms in x(6) coincide for this model and for an 
exactly solvable model suggested by Moshe, Neuberger and Shapiro (Phys. Rev. Lett. 
73, 1497 (1994)), the sub-leading terms ~ 6^ being different. Numerical data con- 
firms our analytical calculation. 

PACS numbers: 02.10.Yn, 71.23.-k, 71.23.Aii 

I. INTRODUCTION 

A. Critical power-law banded random matrices and exactly solvable models 

Recently, there was an increasing interest to unconventional random matrix theories 
(RMTs) that interpolate between the Wigner-Dyson RMT and banded RM with the (almost) 
Poissonian level statistics. One of these models is the power law handed random matrix 



theor}**** for which the variance of the off-diagonal elements takes the form: 

PLBRM: m') = \ J (1) 

2l + (^sin(f|z-j|)) 

Here is the matrix size; a and h are two parameters which control statistical properties 
of PLBRM. The variance (0) is nearly constant inside the band |« — j| < 6, and decreases 
as a power-law function (|V^ijP) ~ ~ fo^' K " j| > b. The case a > 1 corresponds 
to the power-law localization which can be found in certain periodically driven quantum- 
mechanical systems^. If a < 1/2 the spectral statistics of PLBRM approaches that of the 
Wigner-Dyson RMT. The special case a = 1 is relevant for description of critical systems 
with multifractal eigenstate a-*^'^'?'^'^ , in particular for systems at the Anderson localization- 
delocalization transition point. On the other hand, it has been conjectureci^ that the spectral 
statistics of critical PLBRM with large h can be mapped onto the Calogero-Sutherland 
model (CS) |s] at low temperature T ~ 1/6. According to this mapping instead of the 
spectral problem of random matrices one studies the statistics of one-dimensional fermions in 
a parabolic confinement potential interacting with the inverse square potential (x— a;^)"^ (for 
real off-diagonal elements in PLBRM, or the orthogonal ensemble ) or non-interacting (for 
complex off-diagonal elements in PLBRM with identical distribution of real and imaginary 
parts, or the unitary ensemble). 

However, there is an RMT for which the mapping onto the CS model is exact ^. This 
is the model of Moshe, Neuberger and Shapiro (MNS) The probability distribution of 
the Hamiltonian H in MNS is given by Pifi) = / df/ V^ifi), where 

Pj>(H) ocexp(^-Tr7i2- (^)'Tr([f/,7i][?7,7i]t)j ; (2) 

the matrix U is either unitary ( for Ti from the unitary ensemble) or orthogonal (for Ti 
from the orthogonal ensemble), and is the Haar measure. 

The connection between PLBRM and MNS is especially clear in the unitary case^ where 
the unitary matrix U = M diag{e'^''^'} M'^ can be diagonalized by a unitary transformation. 
Then the variances of Vij = [m'^'HM^ . . in MNS at given phases (pi are: 

MNS : (|Vi,f ) = . (3) 

One can easily see that Eq.Q coincides with Eq.((TI) at a = 1 if the phases ipn = '2.T^n/N 
are arranged as an ordered array on a circle. In general, the MNS model can be considered 



as an extension of the PLBRM model for the case of a random arrangement of phases 
distributed over the circle with the joint probability distribution P{{ip}) jl^ : 

^(l<^})~TT ^ ' ^ ■ (4) 

sm 

The averaged value of an observable A{H), which is invariant under the transformation 
H —>■ M^TiM, can be calculated as: 

Here (A) stands for the averaging over the Gaussian random matrix H with entries having 
zero mean value and the variance given by Eq.Q. 

The two-point correlation function, which follows from Eq.(^ after the integration over all 
but two phases, was calculated by Gaudin with the help of the model of free non-interacting 
fermions with a linear spectrun^ii: 

2 

, s^{^,-v,){N/27i). (6) 



7^2(s) = l- ^ 



{2nby 



logfe2'rfc_l) e"^ + 1 



If \s\ ^ b, the correlation function is almost constant 7^.2 (^|s| ^ 6j ^ 1. There is a repulsion 
between phases at a small scale controlled by b: 7^.2 (|s| ^ 6] ~ {s/bY- 



B. Spectral statistics of MNS and PLBRM 

The level statistics of RMT is characterized by the density of states 

N 

P{E) = {j:^iE- en) ) , (7) 

n=l 

and its multi-point correlation functions. For example, the two-level correlation function 
R{u) is defined as: 

flM= <<"'"fj;|7^'''' ; {{ah)) ^ (ai) - (aHi) ■ (8) 
The Fourier transform of R{uj) is known as the spectral form-factor K{t): 



hoo 



K{t) = / e"^' R{uj) dio . (9) 



We rescale time by the mean level spacing 

introducing the dimensionless time t = tA. In the limit of small time the spectral form- 
factor K{t 0) is linked to the other important spectral characteristics called the level 
compressibility^^: 

X = lim( lim K{t)) . (11) 



The meaning of x is the following: Let us take a window of the width 6E, 6E/A = n N, 
in the energy space centered a.t E = and calculate the number n of levels inside the window 
at some realization of disorder. The level number variance is T,2{n) = {{n — n)^)- The level 
compressibility is by definition the limit 

X = Jim ( lim ^] . (12) 

n— >oo \N~*oo on J 

The level compressibility contains an information about the localization transition: x ranges 
from XWD = for the Wigner-Dyson statistics with extended wave functions and a strong 
levels repulsion to xp = ^ in the case of localized wave functions and uncorrelated levels 
with the Poissonian distribution. The intermediate situation 

< Xcrit < 1 

is inherent for the critical regime of multifractal wave functionsi^. 

The exact expression for the level compressibility in the unitary MNS reads^^: 

LL 1 [1 - exp(27r6)] f l/(Anb) , 6 » 1 ; . 

Xmns = T—^, TTTT^ - S r- ^^^) 

Li+ 1 [1 - exp(27r6)] ^ i _ , 6 < 1 ; 

where Li is the polylogarithm functioni^. One can see that Xmns obeys an inequality 
< Xmns < 1, at any finite b. 

Moreover, the level statistics of MNS and of critical PLBRM are asymptotically the same 
in two limits: 6^0 and b —>■ 00. 

If b ^ 1, the theory of critical PLBRM with a = 1 can be rigorously developed by 
mapping^ onto the nonlinear supersymmetric a-modeU^. One can show that the level statis- 
tics of critical PLBRM approaches the Wigner-Dyson statistics^^^. In particular, the level 



compressibility of PLBRM goes to zero and asymptotically coincides with the compressibility 
for MNS 

6 > 1 ^ XPLBR.M\a=l ^ XMNS - + ^(r^) < 1 . (14) 

47r 

This is because the phase repulsion in MNS is strong at large h. The phases ipi^j form an 
approximately equidistant lattice-like structure^. 

In the opposite case 6^1, the phase repulsion in MNS is weak and the phases (pij 
do not form a regular structure. Disorder in the phase arrangements at a small distances 
lv^i~V^j| ^ 1/^ inay become especially important and, therefore, there is no a pnon evident 
correspondence between critical PLBRM and MNS at 6^1. 

Let us consider the N ^ oo limit of Eq.(0) at a = 1. If 6 ^ 1 the off-diagonal matrix 
elements of such a PLBRM are parametrically small compared to the diagonal ones 

a = l, {el) = ^:^{\V,^\^)^h''n^-j). n^-j) = \j-^^- (15) 

We will refer to Eq. (fTH|) as to the almost diagonal critical PLBRMs. The parameter (3 cor- 
responds to the Dyson symmetry classes: I3coe = 1 for the Gaussian orthogonal ensemble, 
and (3gue = 2 for the Gaussian unitary ensemble. 

This model cannot be mapped onto the nonlinear sigma model as the mapping is only 
justified if 6 ^ 1. At 6 ^ 1, the compressibility of PLBRM and MNS are close to the 
Poissonian value yp = L The leading correction of the order of Oib) was derived in Refs. 2, 

n n 

I16I | using an approximation of two interacting levels first suggested in Ref . p . Surprisingly, 
disorder in the arrangement of MNS phases does not influence x ^ind the compressibility 
for PLBRM and MNS are again asymptotically the same: 

Xp-XPLBRM\a=l^Xp-XMNS-V2'Kh + 0{h'^)<^l. (16) 



C. Formulation of the problem 

A natural question arises as to whether the level rigidity of critical PLBRM and MNS 
coincide at an arbitrary 6 ~ 1. The numerical simulations^^ did not exclude such a possi- 
bility. The main result of the paper is that it is not the case: the sub-leading corrections 
of order 0{b^) are different in those two models. To prove this statement we analytically 
calculate the second coefficient of the virial expansior}^ for the level compressibility for the 



critical PLBRM of the unitary symmetry class and compare it with the exact result Eq.()13|) 
for MNS. 

As the analytical calculation is quite involved we undertook an extensive numerical in- 
vestigation of the same problem and found an excellent agreement with the analytical pre- 
diction. 

The paper is organized as follows: we briefly discuss the virial expansion in Section |n] 
and re-derive Eq. ()16|) as the first virial coefficient in Section IIIII The main result of the 
present paper, namely, the second virial coefficient for unitary critical PLBRM, is calculated 
in Section HVl The analytical result is confirmed by the direct numerical simulations which 
are presented in the Section HV (^1 We end the paper with a brief discussion and Conclusions. 



II. THE VIRIAL EXPANSION 



The virial expansion (VE) is a method that allows to study spectral statistics of a disor- 
dered system described by a Gaussian ensemble of the Hermitian N x N {N ^ 1) almost 
diagonal random matrices which have random independent element o^^i^^ : 

{H,,) = 0; {Hi) -> {\H,^,\') . 

It is an expansion in the number of interacting energy levels. Unlike the field-theoretical 
approach, VE starts from the Poissonian statistics and yields a regular expansion in powers 
of the small parameter controlling the ratio of the off-diagonal elements to the diagonal 
ones ~ 6^ <^ 1. The expansion has been represented by the summation of 

diagrams which are generated with the help of the Trotter formula. A rigorous selection 
rule has been established for the diagrams, which allows to account for exact contributions 
of a given number of resonant and non-resonant interacting levels. The method offers a 
controllable way to find an answer to the question when a weak interaction of levels can 
drive the system from localization toward criticality and delocalization. An example of the 
spectral form-factor has been considered in Ref. for a generic dependence of the variance 
(l-ffj^jP) on the difference i — j. It has been shown that a term of the order of b'^~^ is 
governed by the interaction of c energy levels. VE has been applied to DOS in Ref. jl?! . 

VE has been described in detail in Ref. Here, we repeat only its basic definitions and 
final results which will be applied to the model ()15|). VE deals with the following correlation 



function in the time domain: 

k{N, r) = ^(( Tr e"'^"/^ Tr e^^"/^ )) = Kq^N, t) + hk^{N, r) + h'^K^i^N, r) + . . . (17) 

For the constant mean density of states K coincides with K. However, they are different 
if {p{E)) essentially depends on energy E. By analogy with Eq. pi|) one can define the 
quantity 

Xo^limf hm k{N,T)). (18) 

It turns out that at small h there is a simple approximate relationship between x and xo 
(see Appendix of the paper jlfi|): 

The mean density of states for the Gaussian ensemble of almost diagonal RMs with either 
localized or (sparse) fractal eigenstates is close to the Gaussian distribution of the diagonal 
entriesi^: 

(p(i?) ) ^ iV^exp (^-^) ^ A ^ 1^. (20) 

The corrections to this formula start with terms proportional to }? . Thus with an accuracy 
of 0(6^) the unfolding factor in Eq. ()19|l can be taken as ~ 

Each function Ki is governed by the interaction of the i + 1 energy levels. The pertur- 
bative expansion (fT7j) is valid if the limit limjv^oo(-^i) is finite. This can be secured by a 
separation of scales: the level interaction is effectively large at the distances \lo\ < Qi^t = 6 A 
which are parametrically smaller than the mean level spacing A. Otherwise, VE fails and 
one has to take into account an infinite number of the interacting levels. 

III. LEADING CORRECTION TO POISSONIAN LEVEL COMPRESSIBILITY 

The expansion p7p starts with the Poissonian form-factor Kp 

lim Ko = Kp = l 

N^ca 

reflecting a distribution of uncorrelated diagonal matrix elements. The functions are 
given by power series in a large parameter 

x = N\T\b, N = A-^(xN. (21) 



The first correction bKi to the Poissonian spectral statistics in governed by the interaction 
of two energy levels. The series for the function Ki in GOE and GUE reads: 



K, = 2^7tP Y.i-^)'C';\k)n'^\k) x''-'; (22) 

k=l 

.(2) ,M _ (2fc-l)!! 

k\{k-l)\ 



= (23) 



In Eq. ()22|) we introduce the real space sum which depends on the correlation function 
defined in Eq. (tT5|l : 

^ff(fc)-^E'(-^M)' = ^ + o(W)' E'= E + E; (25) 

rrt ^ m rn=—N m=l 

where ( is the Riemann zeta functional. The 1/A^-corrections in yield the dependence 
of Ki on a parameter \T\b = x/N. 
To derive the compressibility xo 

Xo - 1 + &Xo^\ 

we have to put r = after doing the limit a; ^ oo: 

xil^ = ]im(^\\m{K,{x,rb))y (26) 

The small time limit means that we have to neglect all the 1/A^-corrections in Eq.()25|). It 
is achieved if one substitutes TZ^^\k) = limTv^oo n^^\k) for n^i\k): 

m=l 

It is convenient to insert Eq. ()27|) into the series ()22p and to sum over k prior to the sum- 
mation over m: 



Km(k)=Y.(nm)) . (27) 



Here /o,i(' ■ ■) are the modified Bessel functionsi^. The sum over m converges at m ~ a; ^ 1 
therefore it can be converted to the integral /q°° dm. After this integration we find: 

x|,'^|/3=i = -2; xl,'^|^=2 = -vr. (29) 



IV. CORRECTION TO LEVEL COMPRESSIBILITY OF ORDER 



A. The second virial coefficient for the critical PLBRM 

Now we focus on the term of the order 0(6^) in Eq. ()17j) . which is governed by the 
interaction of the three energy levels. In the unitary case we will be considering below the 
expression for K2 readsM: 

r\ 00 

(3 = 2: k2 = ^ c^'\h, h, h)nNih, k2, h) x2ikMks)-2 ^3Q^ 

(3) _ '2hk2h - hk2 - k2h - hh Tjki - 1/2) Tik2 - 1/2) r(A;3 - 1/2) 

r(A;i + ^2 + fcs - 3/2) r(fci + l) r(A;2 + l) r(A;3 + 1) ' ^ ^ 



nNiih}) = I [(^(m) (^H (^(|m - n|) ) 

^ m.n L 



(32) 



The series in r.h.s. of Eq. (l3()|l is three-dimensional. It cannot be reduced to a product of 
one-dimensional series because of the function r~^(A;i + A;2 + ^3 — 3/2) in the coefficient 
()31|) . We will decouple the sums over the indices ^1,2,3 using an integral representatioi>i^ of 
the function T~^{z) 

' - ' f ^^P(^)dt; (33) 



r{z) 2txi Jci t 

the integration contour Ci is shown on FigH Then we change the order of summations over 
ki and integration over t. The real-space summation which is implied in the function TZm 
has to be done at the last step. We assume further that the sums over m and n converge 
at large values of the summation variables and transform the sums to a twofold integral (in 
analogy with the derivation of Xo^^)- This assumption is verified below (see Section riVB|) . 
Using the identities 

,k^{k-l/2) 



U Tik + 1) - ^""VTT-y 



and substituting 00 for N in the limits of the real space integrals over m and n (in analogy 
with Eq. ()27|l ) we arrive at the following expression: 

K2{X, r = 0) = —^^ - Jj dm dn dt exp(t) P\n\P\,n\ [P\m-n\ - 3 Q\m-n\) , (34) 



Im(t)l 



Ci 



C2. 



a=+0 Re(t) 



FIG. 1: The integration contours for the variable t: Ci : { = a > , g] — 00; +oo[} and 
C2 = C+\JC^ where C_ = { 3f?(t) G] -oo;a], = -0}; C7+ = G [a;-oo[, = +0}. 

One can put a = +0. 



where 



= Vl + ^' ^\^\^^^\Jh)^ (36) 

and we have absorbed l/-\/2 into x obtaining 1/2 as a common prefactor. The integrand 
in Eg. ()H4|1 as the function of t has a branching point t = 0. A branch cut may be drawn 
along the negative semi-axis. Since the integrand in Eq. (j!M|l is zero if |t| — 00 at 3?(t) < 
a = +0 and has poles neither in upper nor in lower half-plains, we can transform the 
integration contour Ci to C2 = C+ U C_ (see FiglH) which encloses the branch cut. Fourier 
transforming the functions (|H3HHB|) and using a scaled spatial coordinate 

J = j/x, 

we find the x- independent expression for K2: 



F,(-M,t) = pdJ^^^^^H^, (38) 



V 



F,(M,t) = rij-^^^^vrriH. (39) 

We have introduced an infinitesimal positive constant r], which regularizes the integrals (jHEl- 
inni) at small distances. We will show that K2{t = 0) is finite at r] —>■ +0 thus proving that 
the small distances do not play a role. 

Let us separate out real and imaginary parts of parts of Fi 2(A^,t G C2) accounting for 
a branch-cut of F12 as functions of t along the negative semi-axis: 



F,,2\tec, = Fi^\M) +tsign{^{t)) Fil2\M) , M = ^ , sign($-''^^ | 1, if t G C+ 



-1, if t e C_; 



F2 



°° , C0S(M J') 71 f°° 

dJ = = / I 

J' ^1 - (J')2 2 Jm 

\-> = pdJ'^2!^^;rri7p = f;«(M) + ; (2^(A/) - J,(A/)) ; (41) 



. ..M..., . /^d,J„(,). (40) 

A J/ /i _ ( j/)2 2 7a/ 



(-) 



^ cosjMJ') 

V J' J I - {j'Y 



2 Jo 



M 



dyHo{y) + [Ci (M) - Ci (r^M)] 



^(-) 



(42) 
(43) 



Here, Jq i are the Bessel functions, Hq _i are the Struve functions, Ci is the cosine integral 
functioni^. We will use the property: 

lim [Ci (M) - Ci [r] M) + log(r7)] = . 

One can see that the integral of the real part of (^Fi — 3 over t is zero due to a 
cancellation of the integrals over C+ and C_. Thus, we may keep only the imaginary part 
of - 3F2) in Eq.dnH): 



- 3F2) ^ zsign(^5(t)) [Gi(M) + 62 (M) 



ei(M) = 3 , .2 



.(+) 



(+) 



(+) 



e2(M) = 3{Fl 



F^ 



(+) 



.(+) 



6^(-)^(-)^(+) 



(44) 
(45) 
(46) 



We insert Eq. ljl^ into the expression (jHTjl 

I dMFl[Fi-?,F2) J^dtexp{t)^2 dM (Bi + 62) dtViexp{-t), (47) 



'C2 

and integrate over t obtaining 



(2) 



3^3 

Let us consider the first integral in r.h.s of Eq. 



^hm^ [J^°°(iM (ei(M) + 62(1/)) I . (48) 



-00 



Ti = / dMQi{M) = (49) 
l^^^dM js (IJdy Jo(2/))' (25(M) - MM)) - 2 (|Jdy Jo(y))'| • 

We note that the integral Eq. ()49j) does not contain the regularizer t]. The first term with the 
^-function function can be immediately integrated using dy Jo{y) = 1 [3]. The other 
two terms can be integrated by parts with the help of the standard integrals containing the 
Bessel functions^'=. The result reads: 

J, = ^ fdM (m [Jo {M)f) = Af. (50) 

Finally, we have to calculate the second integral in r.h.s of Eq. 

I, ^ pMe,(M) = 6 (f)'^lim 

^ dAf I I" Ay H„(y) - \og(i,)j ('^^1^ - *(A/) + £ dy My)) + 

/ rM \ roo 

+H_^{M) U dyHo{y)-\og{7])\ j^^ dy Uy) 



(51) 



Unlike the integral Ti, the regularizer of the small distances rj enters the expression for 
the integral X2. We have also introduced the upper limit of the integration over M before 
integrating Eq. ()5ip by parts. At intermediate stages, the boundary terms of the integration 
by parts, which result from the different parts of the integrand in r.h.s. of (j3T|) . diverge in 
the limits — > and ^ 00. However, the diverging contributions exactly cancel out 
at the end so that the final answer for I2 does not depend on rj and is finite in the limit 
Vl — * 00. For example, the coefficient in front of log^(r7) 



lim / dM - m) + r dy My) 



is zero because 

foo coo 

1 



roo POO 

/ dyMy) = 2 dy6{y) 
Jo Jo 



lim I / dM rdyJo(y)\ = - lim l Q ( T dy '^^] \ oc lim ^ = 0. 
The cancellation of log(r7) can be checked in a similar way. We have thus proven that the 



regularization of the Fourier images ()38|39|) does not affect the level compressibility. 

We skip a lengthy intermediate integration by parts and present only the answer for I2: 

X2 = Svr^ / dM (Ho (M) Jq (M)) - — / dM [MHo (M? Jq (M)) . (52) 
JO 4 JO ^ ^ 

To our best knowledge, the integrals of a combination of the Bessel function and the Struve 

function in r.h.s. of Eq. (j52|l are not included in the standard handbooks. We describe their 

calculation in Appendices. Here, we give the results: 

roo 1 

/ dM{Ho{M) Jo{M)) = -; (53) 

JO / 

°dM {mHo {Mf Jo (M)) = ^ - -L^ ; (54) 

^X2 = ^(V3-2), Jl+T,= (2-v^)v^vr^ ^ = (2 - v^) (55) 

We insert formulae ()29l55p into Eq. p9p and obtain the expression for the level compressibility 
of the critical unitary PLBRMs: 



(3 = 2: X = 1 - V2 (vrfe) - 1 (2 - v^) {-nhf 



+ Oih") . (56) 



B. Characteristic spatial scale that governs the compressibility 

Let us estimate the characteristic spatial scale that governs the second virial coefficient 
Xq^\ Firstly, we note that all integrals over t and M converge at \t\char ~ 1 (see Eqs. ljlTjl 
and (UHl)) and Mchar ~ 1 (see Eqs. fl^Uj) and (jSSHSlI))- Returning to the spatial variable 
j = xJ, we may estimate its characteristic scale: 

jchar ~ xJchar ~ x/ Mchar ~ X\J \t\ char / M char ~ X > 1 . 

Therefore, Xo^'' is governed by the large distances m,n, (m — n) ~ a; ^ 1 (see Eq. ()34|1 ) . We 

(2) 

have verified the self-consistency of our calculation scheme for Xo > namely, the replacement 
of the real space sum by the integral is justified. 




FIG. 2: A comparison of the analytical result (|56|) with the direct numerical calculation of x 
b = 0.1: stars mark the result of simulations at the different matrix size; solid line is an interpolation 
that yield xIa^^oo = 0.595 it 0.005; dashed horizontal lines gives the analytical result x ~ 0.605 
(|56jl for b = 0.1; dotted line presents the compressibility of unitary MNS. Inset: the coefficient C2 
of the polynomial fitting Eq. (|57j) for S2 as a function of 1/A^. 

The first virial coefficient Xo^'' is also governed by the large distances of the order of x (see 
Eqs. fj28|) ). We may conclude that the small distances do not contribute to the compressibility. 
One important consequence is that the level compressibility is not sensitive to the periodicity 
of the boundary conditions: If we recalculated Xo ' using the spatially periodic variance 
instead of |T5|) we would again arrive at the same results ()29p55|l . 

C. Numerical test of the results 

A comparison of the analytical result (j56|) with the direct numerical calculation of x is 
presented on Fig|21 The data corresponds to Eq. (0) with 6 = 0.1 for the unitary symmetry, 
Pgue = 2. 



Numerical calculation were done in the range 200 < N < 10 000. The true value of x is 
found from the extrapolation of numerically obtained x{^) to ^ oo. 

The number of realizations ranges from 30 000 for the small matrix size (A^ = 200), to 
200 for the larger one (A^ = 10 000). The level compressibility has been obtained from the 
level number variance Ti2{n) = {{n — n)"^): S2(n) has been calculated in a small energy 
window 6E ^ 0.4 at the band center—. Then, the iV-dependent spectral compressibility 
x{N) has been obtained for each value of as a coefficient of the linear term of a 2nd 
order polynomial fit for S2(n) Q|: 

S2(n, N) = Co{N) + x{N)n + C2{N)n'' . (57) 

The fitting range is 1 < n < 100. To increase the accuracy, we have done both disorder and 
spectral averaging of the data. We have also checked that a small change in the window 
width or in the fitting range practically does not change the results. 

For the critical RMT, it is expected^^ that both x ^-^id C2 have a pronounced (oc 1/A^) 
A^-dependence refiecting the finite size effects in T,2{'n,N). The true value of x is obtained 
at the intersection of the the solid line with the vertical axis at the point 

;^|^^^ = 0.595 ±0.005 (58) 

(see Figj21) which should be compared to the analytical results for the critical PLBRM and 
the MNS exact result given by Eq.(|^ and Eq.lfT^. respectively: 

lim XPLBRMib = 0.1) = 0.6056; lim XMNs{b = 0.1) = 0.6548. (59) 

Given an expectation that the missing term of order ~ 0.001 is negative (as well as the 
term of order b) and thus Eq. (j56|) overestimates x by several parts of 10^^, the numerical 
result is very close to the analytical one and is clearly different to the MNS exact result. 

The true value of C2 is expected to be zero—. Indeed, the numerically obtained coefficient 
C2 goes to zero as A^ increases; see the inset of FigI21 This behavior of C2 confirms good 
quality of our numerics. 



V. CONCLUSIONS AND DISCUSSION 



As we have already mentioned in the Introduction, the level compressibi 
critical PLBRM and of MNS are asymptotically the same both at 6 ^ 1 



with the accuracy up to the terms of order 1/6 and b, respectively. While such a coincidence 
is natural at 6 ^ 1, its origin for 6 -C 1 still remains unclear. 

The main result of the paper, Eg. ()5()j] . is an analytical calculation of the level compress- 
ibility for the critical PLBRM ensemble up to the terms of order 6^ and its comparison with 
the corresponding formula for the MNS model of the unitary symmetry class. 

Our result Eqs. ()56|) shows that the compressibility in MNS is larger compared to PLBRM 

/2 

/? = 2: x|mns-x|plbrm^^(4V3(1 + V2) - (11 + 3^2) ) (Trfo)^ 

in agreement with the numerical simulations for PLBRM, see Fig|21 It is also important that 
the result Eq. ()56p is expressible in a simple algebraic form. The fact that all the intermediate 
sums and integrals are exactly doable in terms of elementary functions is not trivial and may 
indicate that the PLBRM theory is exactly solvable. 

Thus, we conclude that the level compressibility for PLBRM and MNS are not identical 
though very close to each other—. It is tempting to assume that the coincidence is a conse- 
quence of a certain relation between and which holds for both models, so that 
the asymptotic coincidence in spectral statistics for 6^1 automatically leads to that for 
6^1. This scenario can also be related with the existence of a field-theoretical description^li 
which is dual to that of the nonlinear sigma model. 

Based on the leading terms in the b and 1/6 expansion one can guess a possible form of 
a relation between x(6) and which can be then checked using the 6^ and 1/6^ terms. 

From this viewpoint our result Eq. ()56|) is also a very useful step. 

However, maybe the most important conclusion we may draw from the above considera- 
tion is that the virial expansion method^^ is working and helps to obtain solutions to very 
non-trivial problems. 
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Appendix: Integrals containing product of Bessel and Struve functions 



In this Appendix we compute the integrals 

/•OO /"OO 

Inti = / HQ{x)Jo{x)dx , Int2 = / xHn{x)Jo{x)dx. (60) 
Jo Jo 

Integral Inti. Using the integral representation of the Struve function Hq we convert 

Inti to the following form 



2 dt ( r°° r 

Inti = — lim / . ^ ( / dx Jo(x) exp (—a + it)x 

71 o^+O Jo Vl — \J0 L 



2 , /-i dt 
— lim / , — 

71 Jo y/l _ t2 

The inner integral over x is zero^^ at 0<t<l, a = and diverges at t = 1, a = 0. 



J dx Jo{x) sm(t x^ 



0, if < t < 1 
OO, if t = 1. 



Thus we see that Inti is determined by the integration over a small vicinity of the point 
t = 1. The infinitesimal parameter a has been introduced to solve an uncertainty ||0 x cx3|| 
with zero coming from the phase volume at t — > 1: We calculate the integral over x keeping 
the finite a: 

Inti = -- lim f'^Q I I ^ - r , ^ (v/-(^ + ^)| , (61) 

t' = 1 - t, z = -. 

a 

The answer (IK!^ results form (jF)T|l after the substitution 



— (2; + i) = V + I exp (z [tt + arctan(l/z)]) . 
Integral Int2- Let us introduce two auxiliary three-fold integrals 

noting that 

r(l) ^ r(2) _ r^^^ _ 2 _ (2) _ ^(1) _ _ ^(2) _ 2 



/^^ + /^^ = dggJo^(g) = Int^ = - = 2/. 



TT 



(2) 

The idea as how to calculate I2 is very similar to the calculation of Ji: we use the 
propertyi^ 

0, if - l<x-y<l, 



dggJo(g) cos(^q[x - y] 



00, ii X — y = ±1, 



introduce an infinitesimal regularizing parameter 



/f ^ = ^ lim , , , , ^ 

vr^ a^+o J Jo a/1 — a/I — y 

integrate over q keeping the finite a: 



dx 



dy 



3ft 



^ dqqJo{q)exp(^-aq + iq[x -yfj'^ , 



(63) 



r(2) 



lim 



1 dx 



dy 



TT^ 7 Jo a/1 — a/1 — 



3ft 



3/2 



(64) 



^1 + (-a + 2[x - y]^ 

and consider only contribution of two small regions {l — x<^l,y<^l} and {x ^ 1, 1 — y ^ 
1}. After lengthy but rather simple algebra we obtain 

r{2) _ 1 



71 



arriving at the answer (jHH). 
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